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Abstract 

In this paper some critical aspects of the behaviour of breaking lattices subject to 
slow driving forces are briefly reviewed. In particular, fluctuations in the response to 
the variation of external parameters are discussed. 



1 INTRODUCTION 

Antonio Coniglio's contribution to the undertstanding of rupture phenomena in disor- 
dered systems is twofold. On one hand he has contributed indirectly to the field, by fur- 
nishing fundamental knowledge indispensable for dealing with dilute networks, like e.g. 
the exponent for the scaling of the density of cutting bonds in a percolating lattice [Q. 
On the other hand he has also pointed out important properties of the response of these 
systems, such as the multifractality of currents in diluted[|| and two component|| re- 
sistor networks. These studies have continued in the observation, investigation and 
explanation of a large quantity of scaling properties in breaking lattices^, and the 
great amount of work done in this field has shown that scale invariance is a fundamen- 
tal property of rupture and breakdown phenomena of disordered systems Q ||. Indeed, 
disorder plays a fundamental role in generating this kind of behaviour ||. 

Scale invariance has been found to be also a more general characteristic of the dy- 
namical response of these systems to the variation of some external parameter, and 
also a characteristic of the fluctuations of their response. For this reason, the adjec- 
tive "critical" is often used, in analogy with the fluctuations of thermal systems at the 
critical point. In the following we shall mainly focus on this aspect of the breakdown 
phenomena and will consider in particular those situations in which an external pertur- 
bation (driving) is applied slowly with respect to the characteristic relaxation times of 
the system. After briefly reviewing some experimental results (Sec 2.), we shall describe 
how some of the observed features are reproduced by lattice models (Sec 3.), and some 
current points of view on the subject are (Sec 4.). 



2 SOME EXPERIMENTAL FACTS 



duetoMogi§. With the aim of verifying the validity of the Gutenberg Richter law also 
at scales much smaller than those involved in an earthquake, he designed and performed 
some original experiments where a controlled pressure was applied to disks made of a 
mixture of resin and hard grains. He observed that, by increasing the pressure, elastic 
waves were released from some localised region within the sample. By recording the 
maximum wave amplitude of each series Mogi was able to show that the relation between 
the observed amplitudes and their frequency of occurrence was of the algebraic type, 
the one known as the Ishimoto and Iida's law in the seismological field. These findings, 
besides to show in a quantitative way that earthquakes and fractures have common 
features, demonstrated the intrinsic critical nature and the importance of disorder in 
the response of a medium to the external perturbations. 

Laboratory experiments can well reproduce not only the critical features observed 
in distribustion of the energy released by earthquakes, but also those observed in the 
distribution of times. In 1968 Scholz succeeded in reproducing Omori's power law 
describing the number of aftershocks observed in a time t after a main event 0]. By 
analysing the acoustic emission from a fracturing basalt rock he also coonfirmed the 
validity of the Gutenberg- Richter law||. In subsequent investigations Hirata[]l(| showed 
that Omori's law holds also for generic microfracturing processes, at least for times large 
enough after the main event. 

Thanks to the advances in technology and informatics, experiments in this fields 
continue to bring new results. By collecting and analyzing acoustic emission signals 
from concrete-like samples, we have shown JlT[ that power laws not only describes the 
frequency distribution of maximum amplitudes in the aftershock series, but also the 
amplitude distibution of the entire time series, as well as the time lags distribution 
between consecutive meaningful events of acoustic emission. We have found in particular 
that the energy (proportional to the squared amplitude) is distributed according to 
N(E) ~ E~ b with S around 1.3, while for time lags one has P(t) ~ i - ^ with £ approx 1.6 
C3j. In addition, lacking of characteristic scales in the fluctuations of response has been 
brought into evidence by measuring the power spectrum and the autocorrelation of the 
time series. The power spectrum has been found of the l// 7 form, with 7 ~ 0.6jll], fl2|| . 

In more recent experiments Ciliberto and coworkers fi~3|, |l~4|| found S ~ 1.25 for the 
energy scaling and conjectured its universality, whereas £ was found to depend on the 
applied external stress. Acoustic emission recorded during pressurization of spherical 
tanks yielded still more evidences of the critical nature of breakdown phenomena and 
revealed the presence of logarithmic oscillations in the power lawsjlj|. Finally, Maes 
and coworkers |i~6| observed scale invariance in acoustic emission amplitudes, time lags, 
and spatial distance between consecutive events also in a cellular glass, where they 
found C = 1.3, but (5 = 2. 

Another important evidence of the critical response of disordered media concerns 
the roughness of the fracture surfaces |l7]]. After Mandelbrot and coworkers measured 
self-affinity properties in the fracture surface of some steels, similarity in the exponents 
(« 0.8) characterizing many different materials has been pointed out(l^, [l9|]. Further 
experiments have shown later that also 0.5 is observed[|o|. According to some evidences 
lower values seem to characterize slow producing cracks and roughness at small scales, 



3 LATTICE MODELS 



Existence of critical fluctuations in the response of model systems was firstly observed |22[] 
in the fibre bundle model p4| , where an external applied stress is evenly shared by a 
stretched bundle of fibers. The elastic modulus is the same for all the fibers, but each 
one can stand a different, finite, amount of stress. The system is initially unloaded, 
then stress is applied in order to break the weaker fiber of the bundle. The excess stress 
is shared by the other fibres, that therefore become more prone to break; whenver some 
fibre exceeds its own failure threshold it also breaks and stress is redistributed again. 
When no more fibres break stress is increased again, up to break the weakest of the 
survived fibers, and so on until all the fibres of the bundle are broken. Each breaking 
process is carried on at constant stress. By exact calculation Hansen and Hammer were 
able to show that as the applied stress is increased from zero to the global failure value, 
there is a probability P(s) ex s~ a that s fibers break in correspondance of the same 
value of applied stress. The exponent was found to be a = 5/2 and largely independent 
of the statistical distribution of fiber strenghts. Hansen and Hammer also observed 
power law distributions for "avalanches' of broken fuses in the numerical simulation 
of a square resistor network where, in analogy with he fiber bundle model, each fuse 
posseses the same conductance, but can stand different maximum current. When a fuse 
finds itself above its own threshold of rupture, it breaks and the excess current is shared 
by the other fuses. It may happen therefore that some other fuse burn in turn, and so 
on. They found numerically a = 2.7 for this system, very close to 5/2. 

Presence of similar behaviour was later observed in system with vector elasticity by 
numerical simulations of the Born model on a triangular lattice |2f|. In this model sites 
interact via the potential 

Vij = (a - (3) [(iii - itj) ■ ry] + (3 [Ui - Uj] , (I) 

where Ui is the displacement vector of site i from equilibrium, fij is the unit vector be- 
tween the initial equilibrium position of sites i and j, and a and (3 are force constants. In 
the case of Ref. [ p5| avalanches were also triggered by a corrosion mechanisms, accord- 
ing to which the bonds neighbouring a broken bond are weakened. Such a mechanism 
enhances the critical properties of the system[p6|, and the related exponent was found 
to be a — 2.0, the same found in the experiments on the cellular glass Jl6t . 

The power law distributions mentioned above are computed by considering all the 
avalanches occurring during the life of the system. One can also consider what happens 
for a given value of the external solicitation (stress, strain, current) a when averaging 
over many realizations. For the fibre bundle it has been shown that the following scaling 
form holds HI: 

P(s,a) = f(±)s- T (2) 
so 

where 

s « (a c - a)- K , (3) 

and t — 1.5 and k, — I; a c represents the critical value of stress at which the network 
definitively tears. Thus, if one cumulates all the avalanches from the beginning, a = 0, 
to the end of the system life, a = a c , one observes 



with a = t + 1/k = 2.5. Zapperi at al. 31 succeeded to show that within the effective 



medium approximation the fuse network obeys the same scaling of the fibre bundle 
model, and that the ensemble averaged burst size (s) diverges at criticality as 

(s)K(a c -a)^ (5) 

with 7 = 1/2. These results have been confirmed numerically both on square fuse 
networks |3lJ and on elastic networks with central and bond bending potentials |32[ . 

A way for quantifying the vicinity of a system to some critical point is to define a 
branching ratio p for the process as the probability for the breaking of a bond to give 
rise to the breaking of another bond. In analogy with other processes ||| at the critical 
point p = 1. We have evaluated p for the Born potential on a triangular lattice (2(| 
by averaging over many realizations the number of broken bonds at fixed values of the 
external driving (stress or strain) (s). This quantity is related to p by: 

and therefore 

pKS l-(a c -a)-*. (7) 

Numerical simulations show that different types of driving produce different behaviours, 
as expected both from experiments and models Q. By computing p for the Born model 
we have shown in particular that breaking the bonds by driving with strain makes the 
system less critical than breaking them by driving with stress, and that the presence of 
corrosion mechanisms leads closer to criticality. 

It must be stressed that the scale invariance for the breaking of bonds in model 
systems cannot be directly assimilated to that of acoustic emission. In fact, the latter 
carries an energy content that is not included in the broken bonds counting. Self- 
similarity of the energy bursts in lattice models only holds at comparable values of the 
external driving applied |25|, or when special conditions are imposed to the system |34[|. 

Self affine properties of the fracture surface are rather well reproduced by model 
simulations. In particular, numerical results suggest that small self-affinity exponents 
(around 0.5) are related to slowly moving cracks, whilst exponents close to 0.8 are 
associated with fast cracks Some different exponents, like 0.42 have been also 

found for slowly developping cracks in two dimensions \ ^T\ , but also if simulations well 
reproduce generally the experimental observations, there is no a comprehensive theory 
for all these behaviours (for a recent review on dynamic fractures see Ref. |3§| l . 



4 POINTS OF VIEW 

A coherent description of breakdown phenomena and of their critical properties is lack- 
ing. One possibility is to relate them to the Self Organized Criticality (SOC)(27j]. In 
this picture the SOC state corresponds to the final breakdown. Another tempting 
possibility is that of exploiting the analogies with critical phenomena for setting up a 
thermodynamic- like description of the fracturing process. Within this framework Zap- 
peri et al.pH suggested to choose the elastic moduli (or the conductance for the fuse 



in usual first order transitions, would originate in this case from the presence of long 
range interactions; In different situations these have been shown to give rise to frac- 
tal fluctuations in the nucleation process]^. The mean field picture, from which this 
interpretation is derived, does not describe all the features of the breaking process in 
euclidean lattices: pheraps the most relevant shortcoming is that the theory does not 
predict the scaling of the breaking stress distribution with the system size pit. More- 



over, Sornette [B9L 40(1 and da Silveira [Hi, 42 have found that the choice of different 



disorder distributions can select between continuous and discontinuous transition, with 
a tricritical point separating the two behaviours in the phase diagram of the system. 
Finally, continuous breakdown transition can also be observed when constant strain 
instead of constant stress driving is applied to the lattice |4jj|. 

In conclusion we may say that in spite of the observation of critical fluctuations 
in many breakdown phenomena, the physical mechanisms at their origin are still far 
from clear and that it seems necessary to take into account both irreversibility and 
non self-averaging in order to come to a satisfactory explanation of these complicated 
non-equilibrium processes. 
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